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ABSTRACT 

We have found a new evolutionary path to Type la supernovae (SNe la) which has been overlooked 
in previous work. In this scenario, a carbon-oxygen white dwarf (C-l-0 WD) is originated, not from 
an asymptotic giant branch star with a C-l-0 core, but from a red-giant star with a helium core of 
^ 0.8 — 2.0A/q. The helium star, which is formed after the first common envelope evolution, evolves to 
form a C-f O WD of ~ 0.8 — I.IMq with transferring a part of the helium envelope onto the secondary 
main-sequence star. This new evolutionary path, together with the optically thick wind from mass- 
accreting white dwarf, provides a much wider channel to SNe la than previous scenarios. A part of 
the progenitor systems are identified as the luminous supersoft X-ray sources or the recurrent novae 
like U Sco, which are characterized by the accretion of helium-rich matter. The white dwarf accretes 
hydrogen-rich, helium-enhanced matter from a lobe-filling, slightly evolved companion at a critical rate 
and blows excess matter in the wind. The white dwarf grows in mass to the Chandrasekhar mass limit 
and explodes as an SN la. A theoretical estimate indicates that this channel contributes a considerable 
part of the inferred rate of SNe la in our Galaxy, i.e., the rate is about ten times larger than the previous 
theoretical estimates for white dwarfs with slightly evolved companions. 

Subject headings: binaries: close — novae — stars: mass-loss — stars: supernovae — white dwarfs — 
X-rays: stars 



L INTRODUCTION 

Type la supernovae (SNe la) have been widely believed 
to be a thermonuc lear explosion of a m ass-accreting white 
dwarf (WD) (e.g., Nomoto et al. 1997 for a recent review). 
However, the immediate progenitor binary s ystems have 
not been identified yet ( Branch ct al. 1995| ). There ex- 
ist two models discussed frequently as progenitors of SNe 
la: 1) the Chandrasekhar (Ch) mass model, in which a 
mass-accreting carbon-oxygen (C-l-0) WD grows in mass 
up to the Ch mass and explodes as an SN la, and 2) the 
sub-Chandrasekhar (sub-Ch) mass model, in which an ac- 
creted layer of helium atop a C-l-0 WD ignites off-center 
for a WD mass well below the Ch mass. The early time 
spectra of the majority of SNe la arc in excellent agreement 
with the synthetic spectra of the Ch mass models, while 
the spectra of the sub-Ch mass mo dels are too blue to 
be comparable with the o bservations ( Hofiich fc Khokhlov 
1996| ; [Nugent et al. 1997| ). 

For the evolution of accreting WDs toward the Ch mass, 
two scenarios have been proposed: 1) a double degenerate 
(DD) scenario, i.e., merging of double C-fO W Ds with 
a combined rn a ss surpassing t he Ch mass limit (Ibcn & 
Tutukov 1984|; IWebbink 1984| ), and 2) a single degener- 



still debated (e.g., Branch et al. 1995), although theoreti 



ate (SD) scenario, i.e., accretion of hydrogen-ric h matter 
via mass transfer from a bin ary companion (e.g., Nomotc 
1982at [Nomoto et al. 1994|) . The issue of DD vs. SD is 



cal modeling has indicated that the merging of WDs leads 
to the accretion-induced collapse rather than SN la exp lo- 
sion (Saio & Nomoto 1985, 1998; [Segretain et al. 1997|) . 

For the Ch/SD scenario, a new evolutionary model has 
been proposed by Hachisu, Kato, & Nomoto (1996; here- 
after HKN96). HKN96 have shown that if the accretion 
rate exceeds a critical rate, the WD blows a strong wind 
and burns hydrogen steadily at this critical rate and ex- 
pels excess matter in the wind. The WD increases its mass 
up to the Ch mass avoiding formation of a common enve- 
lope. Li & van den Heuvel (1997) have extended HKN96's 
model to a system consisting of a mass-accreting WD and 
a lobe-filling, more massive, main-sequence (MS) or sub- 
giant star (hereafter "WD-I-MS system"), identified with 
luminous supersoft X-ray sources, and found that such a 
system is one of the main progenitors of SNe la as well as 
a system consisting of a WD and a lobe-filling, less mas- 
sive, red-giant (hereafter "WD-I-RG system" ) proposed by 
HKN96. 

Recently, Yungelson & Livio (1998) have reanalyzed the 
models by HKN96 and Li & van den Heuvel based on 
their population synthesis code and concluded that both 
HKN96's WD-^RG and Li & van den Heuvel's WD-I-MS 
systems can account for only (at most) 10% of the in- 
ferred rate of SNe la in our Galaxy. However, Yungelson & 
Livio (1998) overlooked important evolutionary processes 
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both in the WD+MS and WD+RG systems. In this pa- 
per, we first describe an important evolutionary process 
to form the WD+MS sys tem, which has been overlooked 
in previous works (e.g. , |Di Stefano fc Rappaport 1994 ; 



., pi Si 
; |Yungi 



Yungelson et al. 1996; [Yungelson fc Livio 1998). An- 



other evolutionary process leading to the WD+RG system 
is discussed elsewhere (Hachisu, Kato, & Nomoto 1999, 
hereafter HKN99). In §2, we describe the new evolution- 
ary process to form the WD+MS system. Including this 
new evolutionary path, we will show that the secondary 
(slightly evolve d MS star) becomes a helium-rich star like 
in U Sco (e.g., [Williams et al. 198l|), which is transfer- 



ring helium-rich matter onto the primary (WD). We have 
reanalyzed such a helium-rich matter accretion onto the 
WD based on the optically thick wind theory developed 
by Kato & Hachisu (1994). If the secondary MS star has 
a mass of ~ 2 — 3.5Mq, a WD with an initial mass of 
0.8 — I.IM0 grows in mass to the Ch mass and explodes 
as an SN la. We describe the evolution of such a WD+MS 
system in §3. The new parameter region thus obtained 
is much wider than that by Li & van den Heuvel (1997) 
and Yungelson & Livio (1998). Discussions follow in §4, in 
which we have estimated the realization frequency of our 
WD+MS systems that accounts for about a third of the 
inferred rate of SNe la in our Galaxy. 

2. FORMATION OF A NAKED HELIUM CORE AND ITS 
EVOLUTION TOWARD C+O WHITE DWARF 

In this section, we describe the evolutionary path to 
form the WD+MS system where the secondary (MS star) 
has a helium-rich envelope. This important evolutionary 
path shown in Figure ^ has been overlooked in previous 
work. Yungelson & Livio (1998) applied their population 
synthesis code to the WD+MS systems and found that the 
realization frequency of the WD+MS systems is at most 
one tenth of the inferred rate of SNe la in our Galaxy. 
In their population synthesis code, they consider only ini- 
tial systems consisting of a more massive AGB star with 
a C+O core and a less massive main-sequence star. This 
system undergoes a common envelope evolution and finally 
yields a binary system of a mass-accreting C+O WD and 
a lobe-filling MS or sub-giant star. Their code does not 
include another important evolutionary path, in which a 
more massive component fills up its inner critical Roche 
lobe when it develops a helium core of ~ 0.8 — 2.0M© in 
its red-giant phase. 

2.1. Common Envelope Evolution at Red- Giant Phase 
with a Helium Core 

This evolutionary path from stage j4 to F in Figure ^ 
has been first introduced by Hachisu & Kato (1999) to ex- 
plain the helium-rich companion of the recurrent nova U 
Sco. We consider, for example, a close binary at a separa- 
tion a with the primary of mass Mi^i = 7Mq and the sec- 
ondary of M2,i = 2Mq {stage A). When the primary has 
evolved to a red-giant of the radius Ri forming a helium 
core of mass Mi^e (stage B), it fills up its inner critical 
Roche lobe, i.e., i?i — Rl. Here, Rl is the effective radius 
of the inner critical Roche lobe of the primary, which is 
approximated by Eggleton's (1983) formula, 
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0.6g2/3 + ln(l + 5i/3)^ 



(1) 



for the mass ratio q — M1/M2. The long dashed line 
in Figure shows the p rimary radius R^ agai nst the he- 
lium core mass Mi 30 ( Bressan et al. 1993 ). The ra- 
dius of the primary increases with the helium core mass 
from M130 ^ O.2M0 to ~ I.4M0 until helium burning 
ignites at the center of the helium star. If the primary 
fills its inner critical Roche lobe at a certain Miho, i-e., 
Ri{Mi He) = ^1, the separation of the binary, a, is given 

by 



RliMlMc) 

fiq) ■ 



(2) 



Using this relation, we can plot the initial separation a^ 
against Afi^e as shown by the thick solid line in Figure pi 
Then the mass transfer begins. This mass transfer is 
dynamically unstable because the primary red-giant star 
has a convective envelope. The binary undergoes common 
envelope (CE) evolution {stage C), which yields a much 
more compact close binary consisting of a naked helium 
star of Afi3o and a main-sequence star of M2 — 2Mq 
{stage D). Figure shows the separation af^cE and the 
inner critical Roche lobe radius of the secondary R2 t ^g 
after the common envelope evolution. Here, we assume 
the relation 
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(3) 



with the efficie ncy acE = 1-0 for th e common envelope 



evolution (e.g., Iben fc Tutukov 1984; Iben & Livio 1993 



Yungelson & Livio 1998 J 



After the common envelope evolution, the radii of the 
inner critical Roche lobes of the primary and the secondary 
become R\ ~ 0.36a/,cE and i?2 ~ 0.4a/,cE, respectively. 
Since the secondary radius (shown by the dashed line in 
Figure || for the 2M0 ZAMS) should be smaller than its 
inner critical Roche lobe, i.e., i?2 < -R2 f cei the initial sep- 
aration ai should exceed ~ 80 Rq as shown in Figure 0. 
The upper bound of the initial separation is obtained from 
the maximum radius of the IMq star which has formed a 
helium core, i.e., at « 2i?i,max ^ 2 x 300 Rq. Thus, the 
allowable range of the initial separations for our model is 
80 Rq <^ Oi < 600 Rq, so that Rl - 1.4 - 18 Rq and 
i?^ - 1.6 - 2o"i?0 for the case of Figure |. 

2.2. Transfer of Helium Envelope to the Secondary 

After the hydrogen-rich envelope is stripped away and 
hydrogen shell burning vanishes, the naked helium core 
contracts to ignite central helium burning and becomes a 
helium main sequence star of mass Mi, He {stage D). For 
Ml, He ^ 2M0, its C+O core mass is less than l.OTM© 
which is the lo wer mass limit to the non-degenerate carbon 
igniti on (e.g. 
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1988 for a review) . Then the helium star forms a de- 
generate C+O core, whose mass Mq+o grows by helium 
shell burning. When Afc-i-o becomes 0.9 — I.OMq and the 
core becomes strongly degenera te, its helium envelope ex- 
pands to Rl ~ 1.4 — 18 Rq (e.g., Paczynski 1971a; Nomotc 
1982b) to fill its inner critical Roche lobe again {stage E). 
Helium is transferred to the secondary stably on an evo- 
lutionary time scale of tev ~ 10^ yr because the mass 
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ratio is smaller than 0.79 {q = Mi/N'h < 0.79). The re- 
sultant mass transfer (MT) rate is Mi ^ IO^^Mq yr_^ 
which is to o low to form a common envelope (e.g. , |Nec 



et al. 1977; Kippenhahn & Meyer-Hofnieister 197?) 



ter the helium envelope is lost, the primary becomes a 
C+0 WD of A%D '- 0.9 - l.lAfg {stage F). In Figure | 
Mij^MT — MwD is plotted as a function of Mi^e- Here, 
we assume the relation of (in solar mass units of Mq) 



NL 



WD 



0.2 (Mi,Hc - 0.85) + 0.85, 

Ml, He, 



for 0.85 < Mi.Hc <, 
for 0.46 < Ml He < 

(4) 
for the final degenerate C+O WD mass vs. the initial 
helium st ar mass relation ( redu ced from the ev olutionary 
paths by Paczynski 1971a and Nomoto 1982b). Helium 
stars with Mne < 0.85AfQ do not expand to ~ 2 — lOi?©, 
thus burning most of helium to C+O without transferring 
helium to the secondary. For Mi, He ^ 2.0Mq, the C+O 
core has a mass larger than 1.O7M0 before becoming de- 
gener ate, thus igniti ng carbon to form an 0+Ne+Mg core 
(e.g., Nomoto 1984), thus we do not include their mass 
range. For A/iHe <, O.46A/0, helium is not ignited to form 
a C+O core. 

The secondary increases its mass M2 by receiving almost 
pure helium matter of AAfHe ~ 0.1 — O.6M0, which is a 
function of Mi.ho as plotted in Figure H. Then a helium- 
enriched envelope is formed as illustrated in Figure n^ {stage 
F). After the helium mass transfer (MT), the separation 
increases by 10%— 40%, i.e., af^MT '^ (1-1 — l-4)af,cE '^ 
(4 — 70)i?Q. Here, we assume the conservation of the to- 
tal mass and angular momentum during the helium mass 
transfer, which leads to the relation of 



Starting from a close binary of A/wd.o = Mi,c+o wd ^ 
0.8 - I.IMq and AfMs,o = M2,ms ~ 1-7 - S.SA/q with the 
separation of aq = flf.MT ^4 — 80 Rq, we have followed 
a growth of the WD component to examine whether the 
WD reaches 1.38Mq and explodes as an SN la (from stage 
Fto J in Fig. |). 

The initial secondary now has a helium-rich envelope 

{stage F). It evolves to expand and fills its inner critical 

Roche lobe near the end of main-sequence (MS) phase. 

2, Then, it starts mass transfer {stage G). This is a case 

0.85 ^ mass transfer aiter Kippenhahn & Weigert (1967, also 

Paczynski 1971b|) or a cataclysmic-like mass transfer after 
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and we use this to obtain Of mt and then the orbital period 
Po = ^t,MT in Figure |. 

Since the secondary receives ~ 0.1 — O.6M0 helium, its 
hydrogen content in the envelope decreases to X ~ 0.6 
if helium is completely mixed into the central part of the 
star. However, the envelope of the mass-receiving star is 
not convective but radiative so that the helium content 
may be higher in the outer part of the star. 

2.3. Helium,- Enriched Main- Sequence Companion 

We have examined total 5 x 5 = 25 cases, Mi^i = 4, 
5, 6, 7, 9Mq and M2,^ = 1.0, 1.5, 2.0, 2.5, S.OMq, and 
have found the possible progenitors to be in the range of 
Mi,c+o WD ^ 0.8 - I.IM0 and M2,ms ~ 1-7 - 3.5A/0 
with the separation of af^MT ~ 4 — 80 Rq- In these 
cases, the secondary forms a helium-enriched envelope for 
the primary mass of Mi.wD ~ 0.9 — I.IMq (but not for 
Mi^WD ^ 0.8 — O.85M0). We assume, in this paper, that 
the average mass fractions of hydrogen and helium in the 
envelope are X = 0.50 and Y = 0.48 respectively, (assum- 
ing the solar abundance of heavy elements Z = 0.02). For 
the 9M0 + 2.5M0 case, much more helium is transferred, 
while much less helium is transferred for the 6M0 + 2Af0 
case (see Table |l|). 

3. GROWTH OF C+O WHITE DWARFS 



Iben & Tutukov (1984). Since the donor is more mas- 
sive than the accretor (WD component), the separation 
decreases and the inner critical Roche lobe decreases even 
to scrape the envelope mass off the donor star. Thus, 
the mass transfer proceeds on a thermal time scale rather 
than an evolutionary time scale {stage H). The transferred 
matter is helium-rich as observed in the recurrent nova U 
Sco. 

3.1. Optically Thick Winds from Mass- Accreting White 

Dwarfs 

Hachisu, Kato, & Nomoto (1996, HKN96) have shown 
that optically thick winds blow from the white dwarf when 
the mass accretion rate exceeds a critical value {stage 
H). In the present case, the accreted matter to form the 
white dwarf envelope is helium-rich, which is different from 
the solar abundance in HKN96. Assuming X = 0.50, 
Y — 0.48, and Z = 0.02, we have calculated the enve- 
lope models of accreting white dwarfs for various accretion 
rates and white dwarf masses, i.e., Mwd = 0.6, 0.7, 0.8, 
0.9, 1.0, 1.1, 1.2, 1.3, 1.35, and 1.377M0. Optically thick 
winds occur for all these ten cases of Mwd as shown for five 
cases of Mwd = 0.8, 1.0, 1.2, 1.3, and 1.377M0 in Figures 
|5| and ||. Our numerical methods have been described in 
Kato & Hachisu (1994). The envelope solution is uniquely 
determined if the envelope mass Menv is given, where Menv 
is the mass above the base of the hydrogen-burning shell. 
Therefore, the wind mass loss rate Mwind and the nuclear 
burning rate A/nuc are obtained as a function of the en- 
velope mass Menv, i.e., Afwind(Afcnv) and A/„uc(Me„v) as 
shown in Figure pi 

The envelope mass of the white dwarf is determined by 



Menv = M2 - (A/wind + A/„uc). 



(6) 



If the mass transfer rate does not change much in a thermal 
time scale of the WD envelope, the steady-state Afenv = 0, 
i.e., 



M2 = M, 



wind 



Mn 



(7) 



is a good approximation. In such a steady-state approxi- 
mation, the ordinates in Figures || and ^, Af^ind + Mnuo 
are regarded as the mass transfer rate from the secondary 
M2. Thus, the envelope solution is determined from the 
relation in Figure ra for the given mass transfer rate Af2. 
The photospheric radius, temperature, and velocity are 
also obtained from the relations in Figure O. 

Optically thick winds blow when the mass transfer rate 
exceeds the critical rate, which corresponds to the break of 
each solid line in Figure ^. There exists only a static (no 
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wind) envelope solution for the mass transfer rate below 
this break. Its critical accretion rate is approximated as 



Mrr = 1.2 X 10" 



— — 0.40 AIq yr 



(8) 



for X = 0.50, Y = 0.48 and Z = 0.02. If the mass 
accretion rate exceeds the critical rate, i.e., M2 > M, 



For M2 > Mcr, steady hydrogen burning is equivalent 
to the helium accretion at the critical rate of A/cr given by 
equation (||). In this case, wea k helium shell flashes are 
triggered (e.g., Kato et al. 1989| ) and almost all processed 
matter is accumulated on the C+0 WD. Recently, Kato & 
Hachisu (1999) have recalculated the helium wind model 
after the helium shell flashes and estimated the mass accu- 



the styong wind blowe from the white dwarf. The white mulation efficiency with the updated OPAL opacity (ggle 



dwarf! accretes helium almost at the rritiral rate, i e 
-^^nuc ~ -^^cr, and expels the excess matter in the wind 
at a rate of Mwind ~ M2 — Mcr- Here, we assume that the 
white dwarf accretes helium-rich matter from the equator 
via the accretion disk and blows winds to the pole or off 
the equator. 

3.2. Efficiency of Mass-Accretion in Hydrogen Shell 
Burning 

Steady hydrogen shell burning converts hydrogen into 
helium atop the C+O core, which can be regarded as a 
helium matter accretion onto the C+O WD. To estimate 
whether or not the white dwarf mass grows to 1.38Mq, we 
must calculate the mass accumulation efficiency, that is, 
the ratio between the mass accumulated in the WD after 
H/He-burning and the mass transferred from the compan- 
ion. We denote the efficiency by 77H and rjne for hydro- 
gen shell-burning and helium shell-burning, respectively. 
As shown in the previous subsection, the excess matter is 



sias fc Rogers 199(^ ). Here, we adopt their new results in 



a simple analytic form, i.e. 



-5) 



'7Hc 



blown in the wind when the mass transfer rate exceeds the 1999 for more details) 



1, (-5.9<logMHo < 

-0.175(logA'/Hc + 5.35)2 + 1.05, 

(-7.8<logAf < -5.9) 

(10) 
where the helium mass accretion rate, Afuc, is in units 
of A/q yr'-'^. We use this formula for various white dwarf 
masses and accretion rates, altho ugh their results are giv en 
only for the 1.3Mq white dwarf ( |Kato fc Hachisu 19991 ). 

The wind velocit y in helium shell flashes reaches as high 
as ^ 1000 km s~^ ( Kato fc Hachisu 1999 ), which is much 
faster than the orbital velocities of our WD+MS binary 
systems ailorb ~ 300 km s~^. It should be noted that ei- 
ther a Roche lobe overflow or a common envelope does not 
play a role as a mass ejection mechanism because the en- 
velope matter goes away quickly from th e svstem without 
intera cting with the orbital motion (see Kato fc Hachisu 



critical rate, wtiicti leads to rjn = [1^12 — M^ind)/J-VJ-2 < 1, 
for the wind phase. 

During the evolution of mass-accreting white dwarfs, 
the accretion rate becomes lower than Ma in some cases. 
Then the wind stops {stage /in Fig. H). Hydrogen steadily 

burns for A/2 > A/gt ~ 0.5A/cr. Then, we have 77H = 1- 
For M2 < A/st, hydrogen shell-burning becomes unsta- 
ble to trigger weak shell flashes. Once a weak hydrogen 
shell flash occurs, a part of the envelope mass of the white 
dwarf may be lost from the system (e.g., Kovetz fc Prialnili 
1994). In the present study, no mass loss is assumed dur- 
ing the weak hydrogen shell flashes until the mass transfer 
rate becomes lower than A/iow = 1 x 10"'' Af© yr~^, i.e., 
rjn — 1 for A/low <^ M2 < A/st. When the mass transfer 

rate becomes lower than A/iow, however, no mass accu- 
mulation is expect ed from such relatively st rong hydrogen 
shell-flashes (e.g., Kovetz fc Prialnik 1994), i.e., 77H = 
for A/2 < A/low Therefore, we stop calculating the binary 
evolution either when the primary reaches 1.38A/0 {stage 
Jin Fig. |j), i.e., A/i^wd = 1.38A/0 or the mass transfer 
rate becomes lower than Afiow • 

To summarize, we have used the following simplified re- 
lation 

r (A/2 - Afwind)/M2 < 1, (Mcr < M2 < I X IQ-'^Mq 
VH= I 1, (A/low < M2 < Mcr) 

I 0, (M2 < A/low) 

(9) 
for the mass accumulation efficiency of hydrogen shell 
burning. 

3.3. Efficiency of Mass-Accretion in Helium Shell 
Burning 



3.4. Mass Transfer Rate of the Secondary 

We have followed binary evolutions from the initial state 
of (MwD,o, Mms.o, ao) or (Mwd,o, Mms,o, Pq), where Pq is 
the initial orbital period. Here, the subscript naught (0) 
denotes stage F in Figure ^ that is, before the mass trans- 
fer from the secondary starts. The radius and luminosity 
of slightly evolved main-sequence stars are calculated from 
the analytic form given by Tout et al. (1997). The mass 
transfer proceeds on a thermal time scale for the mass ra- 
tio of A/2/A/1 > 0.79. We approximate the mass transfer 
rate as 

A/2 = max ,0 , (11) 

TKH V 4ms / 

where tkh is the Kelvin-Helmholtz timescale (e.g., 
Paczynski 1971b[ ), and (rl = (ilog/?*/dlogM and Cms = 



(ilog/?Ms/c^logM are the mass-radius exponents of the 
inner critical Roche lob e and the main sequence c ompo- 
nent, respectively (e.g., [Hjellming fc Webbink 1987 ). The 
effective radius of the inner critical Roche lobe, R* , is cal- 
culated from equation (|^). 

3.5. Late Binary Evolution toward SN la 
y^ The separation is determined by 

a_ A/1+A/2 ^A/i ^A/2 J .^21 

a Ml + A/2 Ml A/2 J' ^ ' 

We estimate the total mass and angular momentum losses 
by the winds as 



A/ = A/i + A/2 = A/^i„d, 



(13) 
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and 



J 
J 



(Ml 



Af 



(14) 



M1M2 

where ^ is a numeric factor expressing the specific angular 
momentum of the wind, i.e., 



J 2 

-^ = £ ■ a floTb- 
M 



(15) 



For the very fast winds such as Uwind ^ 2a51oib, the wind 
has, on average, the same specific angular momentum as 
that of the WD component, i.e.. 



M2 



M1 + M2 



(16) 



because the wind is too fast to interact with the orbital 
motion (see also HKN99). In equations (12) — ([l5|), we 



must take into account the sign of the mass loss rates, 
i.e., M = Mwind < 0, M2 < 0, J < 0, and so on. 

Figure shows an example of such close binary evolu- 
tions which lead to the SN la explosion. Here, we start 
the calculation when the secondary fills its inner critical 
Roche lobe. The initial parameters are Mwd.o — l-OAf©, 
Mms,o = 2.0Mq, and oq = 9.6i?0 (Pq = 2.0 d). The 
mass transfer begins at a rate as high as M2 = 2.2 x 
10~^AfQ yr~^. The WD burns hydrogen to form a helium 
layer at the critical rate of Mcr = 0.7 x 10~^Mq yr^^ and 
the wind mass loss rate is Mwind = 1.5 x 10~^Mq yr^^. 
Thus a large part of the transferred matter is blown off in 
the wind. 

Since the mass ratio M2/M1 decreases, the mass trans- 
fer rate determined by equation (|ll|) gradually decreases 
below Mcr- The wind stops a,t t — 1.9 x 10"'^ yr. The mass 
transfer rate becomes lower than Mgt ~ 5 x 10~'^Mq yr~^ 
(for MwD = 1.2Mq) at t = 3.8 x 10^ yr and very weak 
shell fiashes may occur. The WD mass gradually grows to 
reach l.SSM© at t = 6.7 x 10^ yr. At this time, the mass 
transfer rate is still as high as M2 = 3.6 x 1O~^M0 yr~^, 
because the mass ratio M2/M1 is still larger than 0.79, 
implying thermally unstable mass transfer. 

This WD-t-MS system may not be observed in X-rays 
during the strong wind phase due to the self-absorption of 
X-rays. However, it is certainly identified as a luminous 
supersoft X-ray source from t = 1.9x10^ yr to 3.8x10^ yr, 
because it is in a steady hydrogen shell burning phase with- 
out a strong wind. Just before the explosion, it may be 
observed as a recurrent nova like U Sco, which indicates a 



helium-rich accretion in quiescence (Hanes 1985|) 



3.6. Outcome of Late Binary Evolution 

Thus we have obtained the final outcome of close bi- 
nary evolutions for various sets of (Mwd,Oi MMS,07ao) or 
(Mwd,o,Mms,Oj-Po)- Figure || depicts the final outcomes 
in the Mms.o — log-Po plane for A/wd.o = I.IM0. Final 
outcome is either 

1) forming a common envelope (denoted by x) be- 
cause the mass transfer rate at the beginning is large 
enough to form a common envelope, i.e., i?i,ph > a '^ 
IORq for M2 ^ 1 X IO^'^Mq yr-i as seen in Figure 



2) triggering an SN la explosion (denoted by ©, Q, or 
0) when Mi,wD = 1.38Mq, or 

3) triggering repeated nova cycles (denoted by A), i.e., 
M2 < Afiow when Afi^wD < 1.38AfQ. 

Among the SN la cases, the wind status at the explosion 
depends on Af2 as follows. 

2a) wind continues at the SN la explosion for Mci- < 
Af2 < 1 X 10"^Mq yr-i (denoted by ®). 

2b) wind stops before the SN la explosion but the mass 
transfer rate is still high enough to keep steady hy- 
drogen shell burning for M^t < M2 < M^ (O)- 

2c) wind stops before the SN la explosion and the mass 
transfer rate decreases to between Afiow < M2 < Afgt 
at the SN la explosion (©). 

The region producing an SN la is bounded by A/ms,o = 
1.8 — 3.2AfQ and Pq = 0.5 — 5 d as shown by the solid line. 

i) The left bound is determined by i?2 f ce = 
i?2(ZAMS) in Figure |, where i?2(ZAMS) is the min- 
imum radius of the secondary at the zero age main- 
sequence (ZAMS). 

ii) The right bound corresponds to the maximum radius 
at the end of main sequence, after which central hy- 
drogen burning vanishes and the star shrinks. 

iii) The lower bound is determined by the decrease in 
the mass transfer rate mainly because the secondary 
mass decreases to reach the mass ratio M2/M1 below 
unity, i.e., Af2 < Miow 

iv) The upper bound is limited by the formation of a 
common envelope. When the mass transfer rate is 
as high as a few times 1O~^M0 yr~^ or more, the 
photosphere of the white dwarf envelope reaches the 
secondary and then swallows it, i.e., i?i,ph <; a (see 
Fig. ^. It may be regarded as the formation of a 
common envelope. The binary will undergo a com- 
mon envelope evolution and will not become an SN 
la. 

The final outcome of the evolutions is also plotted in Fig- 
ures |9|-pl] for other initial white dwarf masses of Mwd.o = 
l.OAfo, 0.9Aro, and O.8M0, respectively. Figure |l| shows 
the regions that lead to an SN la for all the white dwarf 
masses of Mwd.o — 0.75, 0.8, 0.9, and I.IMq (thin solid 
lines) together with AfwD.o — I.OMq (thick sohd fine). 
The region for Mwd,o — 0-7 Mq vanishes. The shrinking 
of the upper bound for smaller Mwd.o is due to larger 
initial mass ratio of Mms,o/Mwd,Oi which enhances the 
mass transfer rate at the beginning of mass transfer {stage 
G), thus resulting in the formation of a common envelope. 
The shrinking of the lower bound can be understood as 
follows: The white dwarf with smaller AfvvD,o needs to ac- 
crete more mass from the slightly evolved main-sequence 
companion. Thus the companion's mass near the stage I 
or J (just before SN la explosion) is smaller after a con- 
siderable part of the hydrogen-rich envelope is transferred 
to supply hydrogen-rich matter to the white dwarf. The 
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thermal time scale of the companion is longer for smaller 
masses, thereby decreasing the mass transfer rate down to 
the nova region. 

4. DISCUSSIONS 

We have estimated the rate of SNe la originating from 
oiir WD+MS systems in our Ca,la,xy hy using equation (1 ) 



of Ibe b & Tutukov (1984), i.e. 



iy = 0.2- Aq- 



Mb 



WT 



dM 

JIM. 



A log A yr 



mary's helium star phase in their scenarios. However, it 
should be noted here that Yungelson & Livio (1998) have 
also obtained a realization frequency of ^-ms ~ 0.001 yr~^ 
under the assumption of no restrictions in their binary 
evolutions, the conditions of which are unlikely. 

A part of our WD+MS systems are i dentified as the lu- 
minous s upcrsoft X-ray sources (SSSs) ( van den Heuvcl ct 
al. 1992). SSSs are characterized by a luminosity of ^ lO'^" 
erg s^^ and a temperature of T ^ 4 x 10^ K {kT ~ 35 
eV), which have been established as a new class of X- 



U' ) ray sources through RO SAT observations (e.g., Kahabka 



& van den Hcuvel 1997 for a recent review). A popula- 



where Ag and A log A denote the appropriate ranges of 
the initial mass ratio and the initial separation, respec- 
tively, and Ma, and Mb are the lower and the upper lim- 
its of the primary mass that leads to SN la explosions, re- 
spectively. For the WD-t-MS progenitors, we assume that 

a,; <; 1 .Snn ffg in order to obtain a, relatively compact con- 



tion synthesis for SSSs has first been done by Rappaport 
Di Stefano, & Smith ( 1994), followed by a m ore complete 
population synthesis (Yungelson et al. 1996). These cal- 
culations predict the total number of the Galactic SSSs of 
~ 1000 and led to the conclusion that the SSS birth rate is 
roughly consistent w ith the observation (e.g., iKahabka & 



dition i of at r.p. after the common envelope evolution. 



van den Heuvel 1997 for a review). Our SN la progenitors 



If the ~ IM0 C+0 WD is a descendant from an AGB 
star, its z ero- age main sequence mass is '^ 7Mq (see, e.g., 
eq.(ll) of Yungelson et al. 1995) a nd the binary separatio n 
is larger than a^ ^ 1300 Rq (e.g., |lben fc Tutukov 1984) ). 
Its separation shrinks to a/,cE ~ 60 Rq after the common 
envelope evolution for the case of acE — 1 and ~ 2Mq 
secondary. Then the orbital period becom es Pq ^ 30 d, 
which is too long to become an SN la (e.g., Li fc van den| 
Heuvel 199?! see also Fig. |). Therefore, the WD-f MS sys- 
tems descending from an AGB star may be rare as pointed 
out by Yungelson & Livio (1998) and may not be a main 
channel to SNe la. 

To obtain the realization frequency of our WD-f MS sys- 
tem descending from a red-giant with a helium core, we 
have followed total ^ 500 evolutions with the different ini- 
tial set of {Mi^i, M2^i,ai) and estimated the appropriate 
range for the initial separation of Alog^ — log Oi, max ^ 
log Oi, mill for total 5 x 5 = 25 cases of (Mi_i,M2,i), each 
for the primary mass of Mio — 4, 5, 6, 7, and 9Mq and 
the secondary mass of A/2,0 = 1, 1-5, 2, 2.5, and 3Mq (see 
Table 0). In Table 0, we omit the case of Mi,o = 4,Mq be- 
cause it never leads to SN la explosions. We find that SN 
la ex plosions occur for the ranges of Af 1 , = 5.5 — 8.5Mr^. 
M-2A h 1.8 - 3.4Af(D, and A log A = 0.5. We thus ob- 



tain the realization frequency of SNe la from the WD-I-MS 
systems as vms = 0.0010 yr~^ for acE = 1, by substi- 
tuting Ag = 3.4/5.5 - 1.8/8.5 = 0.41, Ma = 5.5A//o, 
Mb = 8.5Af0, and A log A = 0.5 into equation (^. For 
comparison, we have obtained a realization frequency of 
SNe la for acE = 0.3. It is still as high as vms ^ 0.0007 
yr~^, which is about one fourth of the inferred rate. Our 
new rate of i^ms — 0.0010 yr^^ is about a third of the in- 
ferred rate of SNe la in our Galaxy and much higher than 
j/Ms = 0.0002 yr~^ (as an upper limit) obtained by Yun- 
gelson & Livio (1998). The reason of their low frequency 
is probably due to the absence of the path through the pri- 



should be observed as an SSS during the steady hydrogen 
shell burning phase without winds, which is about a few 
times 10^ yr as shown in Figure H. Then the number of 
SSSs from our scenario is roughly estimated to be at least 
about ~ 3 X 10^ yr x vms ~ 300, which should be added 
to 1000 by Yungelson et al. (1996). The total number is 
still consistent with observations. 

Our WD-t-MS progenitor model predicts helium- 
enriched matter accretion onto a WD. Strong He II 
AA4686 lines are prominent in the luminous superso ft X- 



ray sources (e.g., Kahabka & van den Hcuvel 1997 for a 



recent review ) as well as in the recurrent novae like U Sco 
i Hanes 1985|; [Johnston fc Kulkarni 1992] ) and V394 CrA 
[Sekiguchi et al. 1989| ). Thus the weakness of the hy- 
drogen emission lines relative to the He II and CNO lines 
is very consistent with the requirement that the accreted 
matter and hence the envelope of the secondary have a 
hydrogen-poor (helium-rich) composition. 

For SNe la, several attempts have been made to detect 
signature of circumstellar matter. There has been no ra- 
dio detections so far. Radio observations of SN 1986G 
have provided the most stringent upper limit to the cir- 



cumstellar density as M/vio = 1 x 10 ^M© yr ^ (Eck et 



al. 1995), where wio means wio = f/10 km s^^. This is 
still 10 — 100 times higher than the density predicted for 
the white dwarf winds, because the WD wind velocity is 
as fast as ~ 1000 km s~^. Further attempts to detect high 
velocity hydrogen signature are encouraged. 
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Table 1 
Initial separations that lead to Type Ia supernova explosions {X = 0.50 and Z — 0.02). 
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''helium mass transferred to the secondary in solar mass units, (minimum, maximum) 
''initial separation of log(ai/i?0), (minimum, maximum) 
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Fig. 1. — Early evolutionary path through the common envelope evolution to the helium matter transfer. 
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Fig. 2. — Long dashed line shows the radius of the primary Ri against its helium core mass Mijjo- The radius increases with the helium 
core mass but it begins to decrease after helium is ignited at the center of the helium core when Mi He ~ 1.4Mq. If the lobe-filling condition 
is satisfied, we obtain the separation a; from equation (2). The thick and thin solid lines show the separations a; (at stage B in Fig. 1) and 
i^/.CE (at stage D) before and after the common envelope evolution, respectively, for the 7Mq + 2Mq pair. The separation shrinks by about 
a factor of ten after the common envelope evolution. Dotted horizontal line indicates the radius of 2Mq star at the zero age main-sequence 
(ZAMS). The helium core mass can grow in mass for the larger separation a;, i.e., for the larger Roche lobe R^ ^. 
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Fig. 3. — Thick solid lines show the mass of the primary (C+O WD) Mi, the mass of the secondary (helium-enriched MS) M2, the mass 
transferred to the secondary AMhoi and the orbital period logPg after the helium mass transfer {stage F) against the helium core mass at 
the beginning of helium mass transfer {stage E). The Dash-dotted line indicates the orbital period at the beginning of helium mass transfer 
{stage E). The separation increases and the orbital period also increases after the helium mass transfer. 
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Fig. 4. — Late evolutionary path to an SN la explosion in our wind model. 
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Fig. 5. — Wind mass loss rate (dashed line) and envelope mass decreasing rate (solid line), i.e., nuclear burning rate plus wind mass loss 
rate, are plotted against the envelope mass for white dwarfs with mass of 0.8, 1.0, 1.2, 1.3 and 1.377M0. The white dwarf mass is attached 
to each line. There exist only static solutions below the break on each solid line while optically thick winds blow above the break. 
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Fig. 6. — Photospheric radius (iJph), temperature (Tpj,), and velocity (tJph) are plotted against the envelope mass decreasing rate (i.e., 
nuclear burning rate plus wind mass loss rate) for five cases of the white dwarf mass, Mwd = 0-8, 1.0, 1.2, 1.3 and 1.377Mq. There exist 
only static solutions below the break on each solid line of the photospheric temperature and radius while optically thick winds blow above 
the break. Photospheric velocity is plotted only for wind solutions. 
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Fig. 7. — Time evolution of a progenitor of SN la. The white dwarf increases its mass to 1.38M0 and explodes as an SN la. Solid lines 
denote the masses of the white dwarf (Mwd) and the main-sequence companion (Mms)- The dashed lines denote, from top to bottom, the 
net mass accretion rate onto the white dwarf (denoted by "WD"), the wind mass loss rate (denoted by "wind"), and the mass decreasing rate 
of the companion (denoted by "MS"), respectively. The dashed line shows the average mass loss rate of ~ 6 X 10~^Mq yr~^ due to helium 
shell flashes after "wind stops." 
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Fig. 8. — Final outcome of close binary evolution in the logPo — Mms.O plane for MwD.o = 1-IMq. Final outcome is either 1) forming a 
common envelope (denoted by x) because the mass transfer rate at the beginning is large enough to form a common envelope, i.e., ijph > lOi?0 

for M2 > 1 X 10~* yr""' as seen in Fig. 6, 2) triggering an SN la explosion (denoted by ©, Qi or ©) when Mi_wD = 1.38Mq or 3) triggering 
repeated nova cycles (denoted A), i.e., M2 < Mi^„ when Mi wD < 1.38Mq. Among the SN la cases, the wind status at the explosion 
depends on M2 as follows. 2a) Wind continues at the SN la explosion for Mcr < M2 < 1 X IO^^Mq yr^-"^ (®). 2b) Wind stops before the SN 
la explosion but the mass transfer rate is still high enough to keep steady hydrogen shell burning for M^t < M2 < Afer (O)- 2c) Wind stops 
before the SN la explosion and the mass transfer rate decreases to between Mi^^ < M2 < Mst at the SN la explosion (0). The region leads 
to SN la explosions is bounded by the solid line. 
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Fig. 10. — Same as in Fig. 8 but for MwD.o = 0.9Mq. 
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Fig. 11. — Same as in Fig. 8 but for MwD.o = O.SMq. 
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Fig. 12. — Regions that lead to SN la explosions arc plotted in the logPo — Mms.o plane for five cases of the initial white dwarf mass, 
i.e., Ms/ifu = 0.75, 0.8, 0.9, and I.IMq (thin solid) together with M^d o = 1-0 (thick solid). The region for M^d o = 0.7Mq vanishes, 
however. 



